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We report the first results of a detailed study of the star formation law 
in a sample of 32 nearby spiral galaxies with well-measured rotation curves, 
HI and II2 (as traced by CO) surface density profiles, and new Ha CCD 



2 I photometry. In this paper we present an atlas of Ha images and radial surface 

brightness profiles and describe a surface- density threshold in the star formation 



law. Prominent breaks in the Ha surface-brightness profiles are identified in 
nearly all of the actively star-forming disks, confirming previous claims of 
star formation thresholds based on lower quality data. We measure the ratio 



^ I of the gas density to the critical density for local gravitational stability at 

^ I the threshold radii. The outer threshold radii observed in Sab-Sdm galaxies 

are in general agreement with those expected from the Toomre Q stability 
I criterion, confirming earlier work, but with a significant variation that appears 

H ' to be weakly correlated with galaxy type. Such a trend could plausibly refiect 

variations in the relative contribution of the stellar disk to the instability of the 
gas disk across this range of galaxy types. Among disks with sub-critical gas 
surface densities, and outside the threshold radius in star-forming disks, the 
number of isolated HII regions increases as the gas surface density approaches 
the critical density. At the thresholds, the gas surface densities span a wide 
range; and the atomic/molecular gas fraction is highest in the disks having 
the lowest total gas surface density. The simple Toomre condition fails to 
account for the active star formation in the inner disks of low-mass spirals such 
as NGC 2403 and M33. An alternative stability criterion based on the shear 



^Data obtained in part at Kitt Peak National Observatory, operated by AURA under contract to the 
National Science Foundation. 
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in the disk provides a better description of these disks but is a less accurate 
indicator of the outer edges of star-forming disks than the Toomrc criterion. 
These results strongly support the view that the formation of gravitationally 
bound interstellar clouds regulates the onset of widespread star formation - at 
least in the outer regions of galactic disks. 

Subject headings: galaxies: ISM - galaxies: evolution - stars: formation - 
galaxies: stellar content 



1. Introduction 

The evolution of galaxies is strongly influenced by how quickly gas is consumed by 
stars. Most simulations assume the star formation rate (SFR) scales as a power law of 
the gas density as originally proposed by Schmidt (1959) or adopt the surface density 
parameterization introduced by Kennicutt (1989). A global correlation between SFR 
and gas surface density {T^sfr oc Hg^s) is obtained empirically when both quantities are 
normalized by the area of the stellar disk. Since galaxies typically have gaseous disks 
that are significantly larger than their stellar disks (e.g. Warmels 1988a; Warmels 1988b; 
Gallagher & Hunter 1984), the gas in the outer disk is not included in the surface density 
measurement (Kennicutt 1998). The empirical star formation recipe implicitly includes 
a surface density threshold for star formation, and this should be incorporated in the 
prescriptions used in numerical simulations. 

The idea that gravitational instability might determine the critical gas density 
for star formation was introduced by Spitzer (1968) and Quirk (1972) soon after local 
stability criteria for differentially rotating disks were developed (Toomre 1964; Goldreich & 
Lynden-Bell 1965). A thin, gas disk is unstable to axisymmetric disturbances where the 
Toomre Q parameter, 

is less than unity. The epicyclic frequency, k, velocity dispersion, a, and surface density, 
H, refer to the gas disk at galactocentric radius R. Widespread star formation is expected 
where the gas surface density exceeds the critical surface density defined as 

^^cr^t = '^Q^' ^ ) 

The parameter ag is fitted to the threshold values of the radially varying quantity 



a{R) = IJ>gas{R)/fJ'crit{R) 



(3) 
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and makes allowances for deviations from the idealized thin disk model such as finite 
scaleheight or the presence of a stellar disk. The velocity dispersion is predicted to remain 
roughly constant in self-regulated regions of disks (Silk 1997) and may have a lower bound 
set by the dissipation of MHD-driven turbulence in the outer disk (Sellwood & Balbus 
1999). If the gas velocity dispersion does not vary much with radius across a spiral galaxy, 
then the critical surface density falls as roughly R~^. This gradient is shallower than the 
decline in the total gas surface density with radius, so the critical density will exceed the 
gas surface density at some threshold radius. Gravitational stability therefore provides an 
appealing explanation for the rather sharp edges of stellar disks (van der Kruit & Searle 
1981a,b; 1982a,b). 

Kennicutt (1989) directly tested this hypothesis using the radial distribution of HIT 
regions to trace the SFR. He adopted a constant gas velocity dispersion of 6 km s~^ and 
found the Schmidt law broke down at radii where the gas surface density was less than 
0.63/icrit (Kennicutt 19890 ). The gas surface density in low surface brightness galaxies 
was subsequently shown to be below the critical density for gravitational instability, and 
gravitational thresholds were used to explain the low level of star formation activity in 
these gas-rich galaxies (van der Hulst etal. 1993; van Zee etal. 1997; van Zee etal. 1996). 
More recently, Ferguson etal. (1998) have questioned the validity of a constant gas velocity 
dispersion. Since measurements of the HI velocity dispersion demand both high angular 
resolution and high brightness sensitivity - requirements not met by most observations, 
the most feasible way forward currently is a comparison of measured threshold radii to the 
thresholds predicted using constant velocity dispersion in a large sample of galaxies. 

Hunter etal. (1998) have recently questioned the utility of the Toomre criterion for 
describing star-formation thresholds. They measure a mean value of aq in irregular galaxies 
that is a factor of two lower than that found by K89 for spiral galaxies. This result implies 
that the gas in the irregular galaxies is less stable than the gas in spiral galaxies. In 
addition to the absence of spiral density waves in irregular galaxies, these disks have less 
rotational shear on average, (i.e. Their rotation curves are closer to solid body.) If the local 
shear rate, rather than the Coriolis force (essentially k), best describes the destruction rate 
of giant clouds, then the stabihty criterion would need to be modified (Elmegreen 1987; 
Elmegreen 1993; Elmegreen 1991). For example, the local shear rate is described by the 
Oort A constant, 

A = -0.5R-. (4) 
Following Elmegreen (1993) and Hunter, Elmegreen, & Baker (1998, hereafter HEB98) 



^ K89 reported a value of aq — 0.67 but used a constant of 3.36 rather than tt in eq. (1). 
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shearing perturbations grow at the rate -kG^/u, and the critical column density for 
significant growth in the presence of shear becomes 



and a A ~ 2.5. This surface density threshold is extremely low in the inner regions of 
galaxies with slowly rising rotation curves but approaches the Toomre stability condition 
where the rotation speed is constant. 

In addition to this degeneracy, many other descriptions of star formation thresholds 
including gas deficiences in the outer disk (Warmels 1988a; Cayatte etal. 1990), the phase 
structure of the interstellar gas (Elmegreen & Parravao 1994), or the radiative feedback 
from young stars have been advanced (Federman ctal. 1979; Skillman 1987). With such a 
broad array of physical processes at work, the K89 analysis should be repeated on a larger, 
more diverse sample of galaxies. The Ha surface brightness profiles presented in this paper 
supercede those in the K89 paper which were based largely on HII region counts. The 
Ha images presented alongside the threshold measurements should resolve the confusion 
which has arisen regarding what is meant by a threshold radius (Ferguson etal. 1998). 
Finally, the number of galaxies with published CO observations has greatly increased in the 
last decade, due in large part to the Five College Radio Astronomy Observatory (FCRAO) 
CO survey (Young etal. 1995). Hence, it is now possible to examine the full range of spiral 
galaxy types using atomic and molecular gas data compiled from the published literature. 

The paper is organized as follows. In §2 we describe our Ha observations and the 
atomic and molecular gas data used in our new analysis. In §3 we identify the star formation 
thresholds, measure gas densities at the thresholds, and evaluate the Toomre stability 
criterion. In §4, we discuss the failures of the model, the shear criterion for stability, the 
infiuence of the stellar disk on the stability criterion, and the implications for the global star 
formation rate. The data presented here will be used to discuss the processes regulating 
the SFR in the high density regime (i.e. where gas surface density exceeds the threshold 
density) in a subsequent paper. 



2. Data 

We selected 32 spiral galaxies divided roughly evenly between types Sab and Sdm 
that had published ^^CO ( J = 1 — 0) and 21-cm radial intensity profiles. Nearly half the 
galaxies are Virgo Cluster members, and 15 of the galaxies are classified as barred galaxies 
in the RCl (de Vaucoulers system). In the K89 paper, radial profiles were estimated for 16 
galaxies from HII region counts. That subsample was dominated by type Sc galaxies. 
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Table 1 lists the type, distance, and disk inclination for each galaxy in our sample. 
The critical density of a disk scales inversely with the assumed distance and is the only 
distance-dependent quantity that enters our analysis. Most of the distances come from 
measurements of Cepheid variables. We correct surface brightness, surface density, and 
rotation speed measurements for disk inclination. Disk inclinations and position angle were 
generally adopted from 21-cm observations. If these values were heavily weighted by HI 
beyond the optical disk, we used the shape of the optical continuum isophotes to derive a 
disk inclination and position angle. The position angles are given in Table 1 and describe 
the major axis of elliptical apertures used in our Ha photometry. 

2.1. Ha Imaging 

We use recombination line radiation from hydrogen to trace the location of young, 
massive stars in the disks. The Ha line was chosen for the imaging because only optical 
CCD cameras offered a wide enough field of view to cover the entire disk in one or two 
pointings per galaxy. The measured Ha luminosities can be directly converted to SFR's 
using the calibration of Kennicutt, Tamblyn, & Congdon (1994) - 

= 1.25 xiy "ergs s- "°"'"°'- ^''> 

This relation was derived using a Salpeter initial mass function, dN{m) / dm = m^^ "^^, for 
stellar masses between m = 0.1 and 100 M0. The stars were evolved along the evolutionary 
tracks of Schaller etal. (1993), and the ionizing flux calculated from Kurucz (1992) model 
atmospheres. The amount of extinction at Ha , Ana — O.TSAy (Miller & Mathews 1972), 
is the largest source of uncertainty in the measured SFR. 

Narrowband continuum and Ha images were obtained at the Steward Observatory 
Bok telescope, KPNO Burrell- Schmidt telescope, and KPNO 36-in telescope. These CCD 
frames were corrected for flxed pattern noise in the standard way using a bias frame, dome 
flat, and sky flat. The 90-in data were obtained with a focal reducer, and scattered light 
produced ghost images which had to be removed in the background subtraction. Foreground 
stars were used to register the images and subtract the continuum emission from each 
narrowband Ha image. The line flux was then corrected for atmospheric extinction and 
the transmission of the filter at Ha relative to the average transmission over the bandpass. 
Data obtained at the Bok telescope and Burrell-Schmidt telescope were flux calibrated 
using observations of standard stars. The flve galaxies observed at the Kitt Peak 36-in 
and NGC 7331, NGC 925, NGC 628, and NGC 4402 were calibrated from the integrated 
fluxes of Kennicutt & Kent (1983) or Kennicutt (1998). The Burrell-Schmidt images do not 
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include [Nil] emission, so we increased their flux by a spectroscopically-determined factor 
ranging from 1.2 to 1.4 to consistently include [Nil] emission. Figure ^compares our Ha + 
[NII]AA6548, 83 fluxes to previously published work. The independent measurements are 
typically consistent to within 15%. 

The radial distribution of Ha emission was measured from surface photometry. The 
optical nucleus of each galaxy defines the center of a set of concentric, elliptical apertures. 
The width of the annuli is a tradeoff between the desire to average over individual Hll 
regions and the need to accurately measure breaks in the profile. A width of about 
10" proved to be a good compromise for most of the galaxies. The profiles changed very 
little if this width was changed by a factor of two. The profiles were corrected to face-on 
orientation. 

To convert the surface brightness profiles in Figure ^ to SFR per unit area, we 
recommend applying an internal extinction correction. Ana = 1.1 mag, derived from 
comparison of the free-free radio emission and the Ha emission from spiral galaxies (e.g. 
van der Hulst etal. 1988). Correction for [Nil] emission in the bandpass reduces the SFR 
by 15% to 30%. Radial gradients in internal extinction are typically not strong compared 
to the scatter in extinction at a given radius (e.g. Zaritsky, Kennicutt, & Huchra 1994, 
Webster & Smith 1983, McCall etal. 1985). An exception is the nuclear region. The central 
star formation rates implied by these profiles are highly uncertain due to the unknown 
amount of dust obscuration, photoionization by an active nuclei, and continuum subtraction 
errors. 

The limiting surface brightness was estimated from the uncertainty in the background 
and the area of the aperture and then verified by adding fake sources to the images. Where 
the covering factor of HII regions was low, photometry could be obtained for individual 
HII complexes using local background measurements. In other regions, variations in 
the mean background level dominate the uncertainty in the photometry. In particular, 
scattered light in the focal reducer camera at the Bok telescope produced ghost images 
that were easily identified by eye but difficult to model and subtract. We used a mean 
background measurement, and then estimated the maximum error using the full range 
of mean background levels on the frame. The errorbars shown on the radial profiles are 
therefore maxima and minima in these regions. 
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2.2. 



Moleculeir Gas Surface Densities 



Only a few galaxies in the sample have been completely mapped in the ^^CO ( J = 1 — 0) 
Hne. FCRAO maps have been published for NGC 5194 (Lord & Young 1990) and NGC 6946 
(Tacconi & Young 1989). The galaxies NGC 5457 (Kenney etal. 1991), NGC 3031 (Sage 
& Westphal 1991), NGC 5236 (Crosthsaite & Turner in prep), and the inner 3.2 kpc of 
NGC 2403 (Thornley & Wilson 1995) have been mapped with the NRAO 12m telescope. 
For the rest of our sample, the CO radial profiles represent measurements made along the 
major axis of each galaxy. These profiles have a resolution of ~ 45". The data are described 
in detail in survey papers by Kenney & Young (1988) and Young etal. (1995). 

These source papers adopt CO-to- H2 conversion factors based on Galactic calibrations. 
We follow this convention but adopt a common scale of 



(Bloemen et al. 1986) for all the galaxies. The beam-averaged brightness temperature, 
Tr, inferred from the measured antenna temperture depends on the source-beam coupling 

efficiency and the forward scattering and spillover efficiency as described, for example, by 
Kenney & Young (1988) and Young etal. (1995). We adopted their calibrations, rjfss = 0.75 
and rjc = 0.73, for the FCRAO survey data. The NGC 5194 profile uses rifss = 0.70 and 
?7c = 0.69 — 0.81 following the model of Lord & Young (1990). The NRAO 12m source-beam 
coupling efficiency derived for MlOl (Kenney etal. 1991) was applied to the M81 data. The 
H2 gas surface density was multiplied by a factor 1.4cos(i) to correct for inclination and 
to include the He mass. 



The HI radial profiles were taken from 21 cm observations listed in Table 1. The 
profiles from the Warmels (1998c) paper were derived from strip scans across each galaxy 
with the Westerbork Synthesis Radio Telescope. The half-power beam width was about 
12" in right ascension by 12" {sin{6))~^ in declination. These scans are essentially major 
axis profiles for many of the galaxies, but the scans are aligned closer to the minor axis of 



NGC 4178, NGC 4321, NGC 4501, NGC 4535, NGC 4548, and NGC 4689. 

The other profiles were obtained by integrating over eUipses and correcting to face-on 
orientation. The typical half-power beam width of the other WRST maps is 30 - 60" (Wevers 
1984; Bosma 1981; Bosma 1977; Bosma etal. 1981; Rots 1975). The WSRT observations 
of NGC 5194 (Tilanus & Allen 1991) and NGC 5236 (Tilanus & Allen 1993) have higher 



N{ H2) = (2.8 X 10^° cm-2) Ico{K{Tr) km s"^) 



(7) 



2.3. 



Atomic Gas Surface Densities 
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resolution. The surface density profile for NGC 6946 comes from a Very Large Array map 
(Tacconi & Young 1986) with 40" resolution. Braun etal. (1994) describe the compilation of 
WRST and VLA obseravations used to construct the high resolution HI map of NGC 4826. 

2.4. Gas Kinematics 

Rotation curves are needed to calculate radial stability profiles for each galaxy. Our 
analysis uses the HI 21-cm rotation curves derived by Wevers (1984) and Guhathakurta 
(1988) and position - velocity diagrams published by Warmels (1988a, 1988b). Velocities 
measured from CO observations were included in the rotation curves for NCC 5194, 
NGC 5457, and NGC 4321 (Kenney & Young 1988; Young etal. 1995; Sofue 1996). These 
CO and HI 21-cm observations lack sufficient spatial resolution to determine the kinematics 
in the inner ~ 45" of the disks. Optical emission-line velocities were used to improve inner 
rotation curves of NGC 4548 and NGC 4639 (Rubin et al. 1999). In general, the dynamics of 
the inner disk are not well constrained by the observations compiled here, and our stability 
analysis will focus on the outer regions of disks. 

Measurements of the HI velocity dispersion are difficult to make in galaxies that 
are inclined enough (relative to our sightline) to accurately model their rotation curves. 
Studies of spiral galaxies find cr = 3 — 10 km s~^ in the outer disk (van der Kruit & 
Shostak 1984; Dickey etal. 1990). A similar range is estimated for dwarf irregular galaxies 
- 7.6 — 11.4 km s~^ (van Zee etal. 1997). The higher velocity dispersions tend to be 
associated with the regions of higher star formation activity (Shostak & van der Kruit 1984; 
Jogee 1999). We adopt cr = 6 km s~^ to facihtate comparison with K89. 

3. Testing the Gravitational Stability Model at the Star Formation Thresholds 

To give the reader an appreciation of the accuracy with which the stability criterion 
can be evaluated, we step through our analysis of NGC 5236. Our Ha atlas is then used to 
further illustrate what is meant by the expression star formation threshold. We describe 
the properties of the gas disks at the thresholds and evaluate the local stability of the disks 
using the single-fluid model. 
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3.1. Case Study: NGC 5236 

Figure |^ shows that the broken spiral pattern of Ha emission from NGC 5236 creates 
a smooth radial surface brightness profile, which is well fitted by an exponential profile 
over a wide annulus. The strongest break in the radial profile defines the threshold 
radius, Rhii = 305". The azimuthally-averaged SFR there is 0.0027 Mq yr^^ kpc^^, but 
the best estimate for the upper limit just 500 pc further out is 100 times lower. The 
Ha surface brightness is also low in the inner disk, where the star-forming regions are 
concentrated along the bar; but the magnitude of the deficit relative to the fit is much 
smaller. In Figure ^, a few HII regions are seen beyond Rhii , but their covering factor is 
much, much lower than in the region R < Rhii- The threshold radius derived from the 
azimuthally-averaged Ha profile is a good approximation to the 'edge' of the star-forming 
disk, but it does overestimate the threshold radius at some position angles (e.g. northeast 
side of NGC 5236) since the shape of the Ha isophotes changes slightly with galactocentric 
radius. 

The middle panel of Figure ^ compares the radial variation in the azimuthally-averaged 
HI surface density (Tilanus & Allen 1993, TA93) to the variation in molecular gas surface 
density along the major axis (Young et al. 1995). More recent CO mapping across NGC 5236 
(Crosthwaite & Turner, private communication) shows the molecular gas is concentrated 
along the bar near the major axis. We find the spatial mismatch (one-dimensional cut vs. 
azimuthal averaging) exaggerates the central depletion of atomic gas relative to molecular 
gas. The azimuthal gas distribution is more uniform further out in this disk. For many 
galaxies in our sample, complete two-dimensional mapping of the gas distribution is not 
available. This example draws attention to the systematic errors introduced by using the 
one-dimensional cuts and suggests they are significant in the inner disk when a gaseous bar 
is present. They do not affect the measurement of the total gas density at the threshold 
radius in NGC 5236. When the Ha surface brightness is plotted against the total gas 
density, a strong break is seen at the outer threshold. This log S vs. log fi relation is well 
described by a power law over most of the disk, but the spatial mismatch between the 
atomic and molecular gas sampling in the inner disk artificially adds some structure to it. 

The critical density for local gravitational instability, as defined by the Toomre Q 
criterion, is also shown in the middle panel of Figure ^ The gas velocity dispersion is 
assumed to be constant, so the small-scale features reflect structure in the rotation curve of 
TA93 through the epicyclic frequency. 

The rotation curves typically represent an average V(R) fitted to the two-dimensional 
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velocity field or to both sides of the major axis. The ratio of the total gas density to the 
critical density, bottom panel of Figure 0, falls steadily with radius beyond the bar region. 
The spatial resolution of the rotation curve is insufficient to describe the critical density at 
radii less than ~ 45" . At the star-formation threshold Rhii , the gas density, 10 M0 pc~^, 
is 70% of the critical density for instability (^(Rhii) = 0.7), as given in eq (2) for o = 
6 km s~^ . In NGC 5236, the dominant error in a (Rhii ) is the extrapolated molecular 
gas density. If the molecular gas disk is truncated abruptly at i? ^ 250" , then the value 
of a (Rhii ) could be as low as 0.25. This scenario is unlikely as it would produce a strong 
discontinuity in the total gas surface density. However, fitting the exponential profile to 
only the less sensitive data {Rco = 150" instead of Rco = 250") yields a fiatter gradient in 
the molecular gas density and overestimates the threshold gas density by ~ 40%. This type 
of systematic error appears to dominate the systematic error in a (Rhii ) for a number of 
galaxies in our sample. 

3.2. Threshold Radii 

Figure |^ shows Ha + [Nil] images and radial surface brightness profiles for the entire 
sample. The threshold radii are marked. Star formation thresholds, as described for 
NGC 5236, are easily identified in 27 of the 32 galaxies in the sample. In most of these 
galaxies the drop in SFR per unit area at the threshold is very sharp, declining by a factor 
of 3 to 400 over a single 10" resolution element. The change in star formation efficiency, i.e. 
SFR per unit gas mass, is similar. The profiles for NGC 2903, NGC 4535, and NGC 4736 
do not show such strong discontinuities, but the sharp change in profile slope in the outer 
disks do seem to refiect some type of threshold in the star formation law. The Ha surface 
brightness in NGC 2841 and NGC 4698 is too low to determine whether the threshold effect 
is present. 

In many disks, HII regions are seen beyond the threshold radius (Ferguson etal. 1998; 
Lelievre & Roy 2000). Their covering factor is much lower than it is within the main 
star-forming disk {R < Rhii)- It seems plausible that the local figas / fi'crit ratio is higher 
in these regions. In NGC 2903 for example, the HII regions beyond the threshold radius 
are clearly associated with density perturbations caused by the spiral arms. The closer the 
mean gas density is to the critical density, the more frequently we expect a given region to 
experience a perturbation sufficient to bring it above the threshold. 
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3.3. Sub-Critical Disks 

Our sample contains seven sub-critical disks where a never reaches the empirically- 
defined threshold value - ag = 0.63 in K89 and aq = 0.69 in a subsequent section of this 
paper. One member of our sample with a sub-critical disk, NGC 2403, shows widespread 
star formation activity, and this failure of the Toomre-criterion will be discussed later. 
Figure ^ shows the instability profiles for the other six galaxies with sub-critical disks. The 
gas density becomes a progressively larger fraction of the critical density (at the radius of 
maximum a (R)) in this order - NGC 4698, NGC 4639, NGC 3031, NGC 2841, NGC 4548, 
NGC 4571. 

To facilitate comparison of their HII region distributions, the six galaxies are displayed 
with the same intensity limits in Figure ^. Starting with the most stable disk in the sample, 
NGC 4698, only a few, small HII regions are detected. The HII regions in NGC 3031 
are confined to the spiral arms whose extent may be set by an outer Linblad Resonance 
(cf. Adler & Westpfahl 1996, but also Westpfahl 1998). The distribution of HII regions 
is much more widespread in NGC 2841 consistent with an instability that approaches the 
critical value. In NGC 4548, the annulus extending from about 0.5 Rhii to Rhii is near 
the instability limit, and a pair of spiral arms is present there. The bar in NGC 4548 
and ram pressure stripping have strongly influenced the gas distribution (VoUmer et al. 
1999) and consequently the star formation. This comparison indicates that - (1) external 
perturbations have an important impact on where star formation occurs in sub-critical 
disks, and (2) HII regions become more prevalent as the instability parameter increases. 

The paucity of star formation in these six galaxies cannot be attributed to low gas 
surface density. The disk of NGC 4698 is a special case in as much as the gas was probably 
removed in a collision with another galaxy (Valluri etal. 1990). The gas surface density in 
3 of the other 5 sub-critical disks exceeds the median threshold density (see next section). 
This result has implications beyond early-type spiral galaxies. Low surface brightness 
galaxies also present a sparse distribution of HII regions and sub-critical gas densities (van 
der Hulst etal. 1993), but the paucity of HII regions can be attributed to the low gas 
surface density (only a few solar masses per square parsec). Our analysis strongly suggests 
that it is the ratio rather than the gas surface density alone that sets the threshold for 
widespread star formation. 
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3.4. The Relative Fraction of Atomic and Molecular Gas 

Using the data described in § we measured the surface density of atomic and 
molecular gas at Rhii in each galaxy. Figure |^ shows the gas surface density ranges 
from 0.7 to 40 M0 pc~^ at the outer threshold. In a few galaxies with prominent gas 
rings, like NGC 4402 and NGC 5055, the apparent edge of the molecular gas disk does 
appear to be associated with the edge of the star-forming disk. However, in half of the 
galaxies in our sample, atomic gas contributes most of the gas surface density at Rhii 
. In eight galaxies, the molecular gas surface density is essentially zero at Rhii for any 
reasonable extrapolation of the CO intensity profile, and threshold radii predicted from the 
molecular gas density alone would grossly underestimate the observed threshold radii. (The 
individual star-forming regions almost certainly contain unresolved molecular clouds, but 
their contribution to the gas surface density is apparently negligible in these regions.) Any 
viable threshold mechanism must explain the wide range in figas (Rhii ) o^nd not depend on 
whether the gas is in atomic or molecular form. The gravitational stability model meets 
both of these criteria. We emphasize that it predicts the formation of large, self-gravitating 
clouds not stars. The good agreement with the observed radii of the star-forming disks, 
however, strongly suggests that star formation will happen once the clouds form. 

At the median threshold gas density the atomic gas fraction ranges from 20% to 100%. 
Yet at densities twice as high, the ISM is almost entirely molecular. The gas is almost 
entirely atomic at surface densities that are half the mean. Much of the gas mass is in 
molecular form in disks with high gas surface densities, including the sub-critical disk of 
NGC 2841. The high atomic gas fractions in the other sub-critical disks are consistent with 
their average to low gas density. This division in atomic gas fraction is expected if the 
column of HI required to shield the molecular gas from radiation (that would dissociate 
it) is fairly uniform among disk galaxies. The implied threshold column for molecular gas 
formation ranges from 5 to 15 Mq pc~^. 

It is interesting to examine whether the disks that are mostly molecular at Rhii become 
dominated by HI where the gas surface density drops into the 5 to 15 Mq pc~^ range. 
Generally, we find that no CO emission is detected at this radius, so the data for the outer 
disks are consistent with the same H2 /HI self-shielding column density. The only disk 
that clearly has a low atomic gas fraction at low total gas column is NGC 4579. 
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3.5. The Instability Parameter at the Outer Threshold 

Critical surface densities for instability were computed from the rotation curves and 
compared to the gas densities described in the previous section. The median value of the 
instability parameter at the outer thresholds, aq = 0.69, was found to be consistent with 
the K89 value. To reveal the physical limitations of this model, we examine the stability of 
the disks at the threshold radii in greater detail than previous work. Three types of second 
order effects are examined: (1) Does the scatter in the a (Rhii ) values reflect the accuracy 
of the model or the magnitude of systematic errors?, (2) Can the sharpness of the threshold 
be predicted from the stability properties of the disk?, and (3) Is the stability of the gas 
disk influenced by the stellar disk? 



3.5.1. Disk Asymmetry 

The epicyclic frequency and, consequently, the critical density have significant errors 
where the rotation curve is poorly determined as in highly inclined disks, disks with 
kinematic disturbances near Rhu , and disks with offset dynamical and photometric 
centers. Within our sample, the projection effects are largest for NGC 628 where a 1° 
uncertainty in the disk inclination implies a 20% uncertainty in ^crit ■ The rotation curves 
of NGC 5457 and NGC 5194 are not well defined at Rhii due to tidal disturbances, so 
the threshold instability parameter is not defined in these two galaxies. The dynamical 
center of NGC 925 is not coincident with the photometric center (in Ha , continuum, or 
HI) - a situation common in barred Magellanic irregulars (Pisano et al. 2000 and references 
therein). In general though, uncertainty about ^cru is small compared to the error in the 
Ugas measurement. 

The largest systematic errors in the measurement of Ugas are introduced by the 
asymmetry of the star-forming disks - i.e. the lack of axial symmetry or roundness. To 
quantify the impact on the threshold measurements, we measured a threshold radius and 
stability profile along the major axis of each galaxy. Comparison to the azimuthally-averaged 
quantities reveals the magnitude of the azimuthal variations in the instability parameter. 
The largest variations are illustrated in Figure |^. We find the following correlations between 
the instability parameter and the distribution of star formation in these disks. 

1. Recent star formation activity in NGC 4254 is skewed toward the north-eastern half 
of the disk. The threshold radius on the western side is only 60% of its average value. 
The instability parameter along the eastern side of the major axis is also higher than 
it is on the western side. The threshold radius and stability parameter derived from 
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the azimuthally averaged profile are heavily weighted by the eastern side of the disk. 

2. In NGC 4321, the average threshold radius Rhii provides a poor description of the 
irregular outer, Ha isophotes. At a position angle of 90° for example, the HII regions 
extend to a radius of O.QRhii to the west but only 0.6Rhii to the east. Inspection of 
Figure |^ shows the two-arm spiral pattern also loses its definition at these radii. The 
surface density is ~ 0.7 ^crit at the threshold along either direction of this cut because 
the gas density declines much more steeply with radius on the east side of the galaxy. 
In contrast, when the gas densities from each side of this disk are averaged together, 
the value of the instability parameter at Rhii is underestimated by a factor of two! 
Better spatial coverage drives the threshold stability parameter toward ^(Rhii) = 0.7. 

3. The star formation activity in NGC 4654 extends furthest to the south-east, but the 
HII regions are faint and sparse here. The Ha surface brightness is much higher on 
the western side of NGC 4654 where the star formation threshold is much sharper. 
The gas surface density is sub-critical over much of the eastern side of the galaxy and 
super-critical on the western side - consistent with the paucity of HII regions to the 
east. On the eastern side, the instability parameter decreases more slowly with radius, 
however, explaining the larger radial extent of the HII regions. The sharp gradient in 
the instability parameter on the western side is likely not resolved, so the value of the 
instability parameter, a = 1.6, at the western threshold is overestimated. 

4. In NGC 7331, the star formation activity is skewed toward the southern half of the 
disk. The azimuthally-averaged Ha profile is weighted toward this side of the disk, so 
the southern edge of the star-forming disk defines Rhii • Along the southern side of the 
major axis, the gas density is 0.8 ficrit at Rhii • At this radius, a = 0.4 on the northern 
side of the disk. However, at the smaller radius defined by the northern edge to the 
star-forming disk, the gas density is O.Gficrit ■ In this case the azimuthally-averaged 
value of the instability parameter, a(RHii) = 0.4, clearly underestimates the true 
value at the edge of the star-forming disk. 

We have demonstrated that using azimuthally-averaged gas densities and SFR's 
can lead to errors in a (Rhii ) as large as a factor of two when the disk is highly 
non-axisymmetric. The theoretical basis for the Toomre stability criterion breaks down in 
these irregular systems, so it is noteworthy that a local Q criterion appears to be so useful in 
describing the distribution of star formation in these highly disturbed disks. The examples 
presented demonstrate that the distribution of star-forming regions clearly becomes much 
more sparse, similar to the subcritical disks, in regions where the instability parameter is 
low. The threshold values of the instability parameter also show a tendency to converge 
when better spatial coverage is used. 
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3.5.2. Star Formation Beyond the Threshold Radius 

Gravitational threshold models only predict a sharp edge to the star-forming disk if the 
gas disk is perfectly smooth. The structure of the interstellar medium in real galaxies must 
modify this picture (Kennicutt 1989). Non-axisymmetric modes and local perturbations 
generate structure on large and small scales, respectively. Such regions may allow cloud 
formation (and subsequently star formation) even where the azimuthally-averaged gas 
density is sub-critical. For example, in § p.3| , the spiral density waves in the subcritical 
disks of NGC 3031 and NGC 4548 are supported where the gas surface density exceeds 
0.3/icrit rather than 0.69/icrit • The lower threshold is consistent with models of swing 
amplification. The more isolated HII regions in the subcritical disks are generated by 
a different mechanism. We suggest that random density perturbations generate some 
super-critical regions in disks that have surface densities close to the critical density 
for instability. Support for this viewpoint comes from the sub-critical disks. Isolated 
HII regions were found most frequently when the gas surface density was only slightly 
sub-critical, and no HII regions were found where the gas surface density was less than 10% 
of the critical density. As the average gas surface density approaches the critical density, 
the amount of star formation activity in isolated, star-forming regions clearly increases. 
Could a combination of local perturbations and spiral arms also explain the existence of 
HII regions at radii beyond the threshold radius? If so, star formation activity should be 
more prevalent where the ^gas I y^crit ratio remains high (i.e. greater than 0.3) beyond Rhii • 

To test this hypothesis, we need to minimize the systematic errors caused by disk 
asymmetry at Rhii , so we compare the 11 galaxies with the smallest azimuthal variations 
in threshold radius. Figure | shows the radial stability profiles for ten members of our 
sample; the profile for NGC 5236 was shown previously. The spread in the instability 
parameter is only a factor of two at the threshold radius where the profiles converge. 
Beyond the threshold radius, the shallowest gradient is found in the disk of NGC 2403 
which remains marginally stable. Not surprisingly, many HII regions are seen beyond Rhii 
in NGC 2403. The gas surface density in NGC 6946 also remains close to the critical 
density, i.e. [igas > O.ificrit, out to 2Rhii ; and a large number of HII regions are seen 
beyond Rhii (cf. Ferguson etal. 1998). In contrast, the steepest gradients in the instability 
parameter, 

I '^^f'sas/^'^"'^ I > 0.1 kpC-\ (9) 

dR 

are found in NGC 4402, NGC 4713, and NGC 5236. None of these galaxies show much 
Ha emission beyond the threshold radius. The instability parameter is also low outside the 
threshold radius in NGC 2903, but the star formation is confined to the spiral arms, locally. 
The likelihood of finding isolated HII regions outside Rhii therefore appears to scale with 
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the value of the instabihty parameter. 

It is interesting to examine why, beyond Rhii , a{R) remains high in some disks and 
drops sharply in others. Figure ^ shows the gas surface density profiles for the same ten 
galaxies. The HI disks of the Virgo cluster galaxies NGC 4394, NGC 4402, and NGC 4689 
appear to have been stripped at radii R > Rhii- When the molecular gas density is included 
the difference in total gas surface density relative to the field galaxies is reduced, but the 
Virgo members still tend to have lower values of of the instability parameter outside Rhii 
. Inspection of Figure ^ shows that the Virgo galaxies show HII regions beyond Rhii less 
frequently than the field sample. 



3.5.3. Variations in the Threshold Stability Parameter with Galaxy Type 

Figure shows our best estimates of a (Rhii ) at the outer edges of 26 disks.0 The 
significance of the weak trend with galaxy type is difficult to assess quantitatively. The 
error bars represent limiting values and are set by the magnitude of systematic errors 
associated with one-dimensional gas profiles in non-axisymmetric disks, tidal disturbances, 
the extrapolation of the CO intensity profile to larger radii, and disk inclination. The 
dominant error terms differ in the early-type and late-type disks, so we discuss their a (Rhii 
) values separately. 



The earlier type disks (Sa to Sb) in Figure |T0] tend to have a (Rhii ) values lower 
than the sample median. The discrepancy is not significant for NGC 2903, NGC 4321, and 
NGC 7731. The range of a values shown for each of them is consistent with the sample 
median, 0.69, because corrections for non-axisymmetric structures were shown to drive a 
in this direction (S p.5.1|) . The a (Rhii ) values for the early-type galaxies NGC 4394, 



NGC 4402, NGC 4579, and NGC 4689 range from 0.5 to 0.6 and depend on extrapolations 
of the CO intensity profile. In all of the galaxies where molecular gas dominates the gas 
density at the threshold, a is less than or equal to 0.69; and most of these galaxies are type 
Sab-Sbc. It is unlikely that the gas surface density is systematically underestimated in 
these disks because the two largest biases in interpretation of the CO measurements work 
to increase the molecular gas density. First, since the CO emission is not well-resolved, 
particularly in NGC 4394, NGC 4402, and NGC 4689, the tendency is to overestimate 
scale lengths and artificially raise the molecular gas surface density at the threshold radii. 
Second, the CO/H2 conversion factor is expected to be lower in disks with low metallicity. 



The six galaxies with sub-critical disks are omitted. NGC 5194 and NGC 5457 are not shown since the 
rotation curves are not well-defined at Rhii ■ 
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If the early-type spirals are more evolved than the late-type spirals, i.e. have consumed 
a larger fraction of their original gas content, then they should have higher metallicity 
on average. Hence, we have no reason to expect a low CO/H2 conversion factor in the 
early-type galaxies. It is plausible then that the fJ^gas/ ^J'crit(Rmi) ratios are somewhat lower 
than the mean in these early-type disks. 

In contrast, the threshold values of the instability parameter for the Hl-dominated 
edges scatter to both sides of 0.69. Lower mean metallicities in the later-type disks would 



introduce an artificial negative slope in Figure |T0|, which is not seen. We find no trend 
toward higher values of the instability parameter in late- type galaxies. In NGC 628, the 
estimated threshold value of the instability parameter, q;(Rhii) = 0.86, is consistent with 
the mean, aq = 0.69, given the uncertainty in the rotation speed described in § ^.5.1| . In 
NGC 925, the instability parameter is overestimated on the western side of the galaxy when 
the offset between the photometric and dynamical centers of this galaxy are neglected. 
(The kinematic radius is smaller than the Ha radius, so k is higher than we estimated 
assuming common photometric/dynamic centers.) A detailed analysis using (jg = 10 km s~^ 
found a is 0.3 at the threshold (Pisano etal. 2000) which corresponds to a = 0.5 for our 
assumed value of cr = 6 km s~^ . Evaluation of the instability is particularly difficult in 
these late-type spirals, owing to the non-axisymmetric structures. We find no trend toward 
lower a (Rmi ) in the low angular velocity galaxies in our sample, all of type Scd to the 
Sdm, as one might have expected based on the low a values reported in dwarf irregular 
galaxies (HEB98). 



3.5.4- Conclusions 

The asymmetries of many disks in the sample compromise the accuracy of the stability 
analysis. We have demonstrated, however, that when disks are examined along position 
angle cuts that eliminate this uncertainty, the convergence of the stability profiles at the 
threshold radii is tighter. Outside the threshold radius, the disks with gas densities well 
below the critical density do not show many HII regions while those with prominent 
populations of HII regions tend to be near the threshold surface density. The weak trend 
toward lower values of ag in earlier type spiral galaxies is consistent with these disks having 
a larger stellar mass fraction (e.g. Jog & Solomon 1984). 
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3.6. Inner Disk Thresholds 

We also examined whether simple gravitational stability arguments might be as 
successful describing star formation thresholds in the inner regions of disks as we have 
shown them to be at the outer edge of star forming disks. 

The sharpest star formation thresholds in the inner disk are associated with rings 
of intense star formation activity - e.g. see NGC 4569, NGC 4579, NGC 4639, and 
NGC 4736 in Figure 0. Does the instability parameter, a(r), increase across these rings? 
The beamwidth of the HI and CO observations is inadequate to evaluate the instability 
parameter across the rings in NGC 4569 and NGC 4639. The major axis instability profile 
across NGC 4579 reaches a = 0.6 on the western side of the ring, but the gas density along 
the eastern major axis is much lower than the critical density. The ring of HII regions in 
NGC 4579 is brighter on the west side than the east side, so star formation in the ring 
appears to be correlated with higher values of the instability parameter, i.e. lower Q values. 
In NGC 4736, the ring of HII regions at i? ~ 45" coincides with an inner Lindblad resonance 
(Gerin et al. 1991; Gu etal. 1996; Wong & Blitz 2000). The spatial resolution of the gas 
data does not resolve the 680 pc wide ring, but the lower limits on the gas surface density 
in the ring are sufficient to raise the instability parameter to values of 0.6 to 1.0. near the 
inner edge of the ring {R ^ 34"). Higher resolution data indicate Q < 1 in localized regions 
of the ring and Q ~ 3 averaged across the ring (Wong & Blitz 2000). 

A few galaxies without rings have gas surface densities less than 0.69/icrit in their 
inner disks. We show the a (R) profiles for the three with the best spatial resolution in 
Figure 0. We expect to see little Ha emission at radii less than the inner threshold radii in 
these galaxies. In Figure ^, the disks of NGC 628 and NGC 4535 show little star formation 
inside the predicted threshold radii of 70" and 50" , respectively. Star formation in the 
central region of NGC 4178 is confined to the bar, where the gas density is likely higher 
than the azimuthal average, so the star-forming regions may not be sub-critical. The inner 
thresholds in these three galaxies are therefore consistent with the gravitational threshold 
model. Only the star formation in the central region of NGC 2403 is in direct confiict with 
the gravitational threshold model. 

The inner regions of many spiral galaxies have a deficit of Ha emission, relative to 
a fitted exponential disk - e.g. NGC 628, NGC 925, NGC 2403, NGC 3031, NGC 4321, 
NGC 4394, NGC 4535, and NGC 5236. Aside from the cases already discussed and 
NGC 4394, these breaks in the Ha surface brightness profile are not large enough to be 
called star formation thresholds. This property of the Ha surface brightness profiles was 
noticed by Hodge & Kennicutt (1983), but the physical cause has not been established. 
The region of low SFR is clearly associated with a bar in NGC 5236, NGC 4394, and 
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NGC 925. The inner threshold in NGC 3031 coincides with the inner Linblad resonance. 
Many galaxies with central Ha deficits have sub-critical gas surface densities in the inner 
disk, and those with supercritical gas densities have bars or resonances. The only galaxy in 
this list with a supercritical gas density in the center and no identified bar or resonance 
is NGC 4321. Apparently, bars and resonanaces can inhibit star formation in the central 
regions of disks that are unstable to axisymmetric perturbations. 

4. Discussion 

Our results support previous studies which have claimed that the outer edges of spiral 
disks are well described by the gravitational stability of a single fluid, isothermal gas disk 
(Kennicutt 1989). We have shown that this criterion is often a good description of star 
formation thresholds in the inner disk as well. The aim of this section is to gain further 
insight into the underlying physical processes that regulate disk star formation. We examine 
the galaxies where the gravitational instabilty model fails and evaluate the utility of a 
quahtatively different threshold condition. A clear prediction of the gravitational threshold 
scenario is that the stellar disks should affect aq. We examine whether first-order accounting 
of the stellar disk changes the scatter in the distribution of a (Rhii ) measurements. Finally, 
the impact of these star formation formation thresholds on the global star formation rate is 
summarized. 

4.1. Shear Models 

4.1.1. NGC 2403 and M 33 

The gas density in the inner 240" of NGC 2403 is below the critical value and less 
than half the critical density within 150" of the center of the disk. Yet the inner half of 
the disk is covered with HIT regions. No other galaxy in our sample presents widespread 
star formation in regions where the gas disk is predicted to be stable. This contradiction, 
and the very similar situation in M 33, were noticed by Kennicutt (1989). The most 
obvious explanations have been dismissed by Thornley & Wilson (1995). In particular, if 
the H2/CO conversion factor were as much as five times higher than the Milky Way value, 
the implied increase in the molecular gas density would bring the total gas density up to 
the critical density in the inner disk without affecting the good agreement at the outer 
threshold. Comparison of virial and molecular masses for clouds in M 33, however, indicate 
the conversion factor is similar to the Galactic one (Wilson 1995). The metallicity and 
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rotation speed of M 33 are similar to NGC 2403, so Thornley & Wilson find no reason to 
suspect a higher surface density of molecular gas in either galaxy. 

The star-forming disks in many dwarf irregular galaxies also appear to be sub-critical 
when the Toomre Q parameter is used to estimate the critical density, and it has been 
suggested that the shear criterion, Eqn. 5, describes the critical density better than the 
Toomre Q criterion (Hunter & Plummer 1996; HEB98). The disks of NGC 2403 and M 33 
share an important dynamical property with dwarf irregular galaxies. The shear rate is low 
in the inner disk where the rotation speed is rising slowly with radius. The disk stability 
analysis for the irregulars, while quite interesting, is not yet compelling. Large systematic 
uncertainties arise from the highly irregular shape of the disks, greater uncertainty about 
the CO/H2 conversion factor at very low metallicities, and the assumption that the 
molecular/atomic gas fraction is constant with radius. 

Figure |1^ compares the shear criterion for the critical density to the Toomre critical 
density across NGC 2403 and M 33. The gas surface densities for M33 data come from 
Newton (1980) and Young etal. (1995), and the threshold radius for star formation is 29'. In 
both spiral galaxies, the gas surface density exceeds /i^^J.^^ but falls below the Toomre critical 
density. The presence of widespread star formation can be interpreted as evidence that 
cloud formation is limited by the time available for cloud growth (in the presence of shear) 
rather than gravitational instability. It remains unclear whether this explanation is correct. 
One disk in our sample, e.g. NGC 4535, has an inner threshold which is well- described by 
the Q criterion but not the shear criterion. Also, models of self-regulated star formation 
predict Q ~ 1 throughout the disk (Silk 1997). Departure from this equilibrium condition 
may indicate the disk is in a transient state. The unusually high infall rate of neutral 
hydrogen onto the NGC 2403 disk, and perhaps M 33 as well, support this view (Schaap, 
Sancisi, & Swaters 2000; Sancisi etal. 2000). 



4- 1-2. Shear Criterion at the Outer Threshold 

The galaxies in our sample generally have fiat rotation curves at the threshold radius, 
so the critical density based on the Toomre Q criterion reduces to 

fJ^crit = OLQ -— . (10) 

ttG 

The shear criterion for the critical density has the same radial dependence where the 
rotation speed is constant, 

fJ'critA = «A ■ (11) 
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Since the fitted value of ag = 0.69, we need = 2.0 to obtain ij,gas = jJ^crUA R-hii for 
tlie disks with constant rotation speed at Rhii . This normahzation factor is very close to 
that predicted by Elmegreen (Elmegreen 1993; Hunter, Elmegreen, & Baker 1998). If the 
rotation speed is constant at the outer threshold, then no distinction can be made between 
these two criteria for the outer edge of a disk. 

The galaxies in our sample with the largest gradients in rotation speed at Rhii are 
NGC 4254, NGC 4394, NGC 4402, NGC 4535, and NGC 4569. The gradient in NGC 4394 
is highly uncertain due to the low resolution of the HI position - velocity diagram for 
NGC 4394. Figure |1^ compares the Toomre (top panel) and shear (bottom panel) 
representations of the critical density for the other four galaxies. The scatter among the 
values of /ig^s //^crit is actually smaller for the shear model than the Toomre Q model. 
However, when the ua term in the shear criterion is normalized using the galaxies with flat 



rotation curves at Rhii , as represented by NGC 2403 and M 33 in Figure |13| the scatter in 
the /igas / ^^cr^t (Rhii ) valucs excecds that obtained with the Toomre Q criterion. The shear 
criterion does not appear to be as robust an indicator of the edges of star-forming disks 
in spiral galaxies as the Toomre criterion. The shear model was also unable to adequately 
describe the edges of the disks in dwarf irregular galaxies (Hunter etal. 1998). 



4.2. The Influence of Stellar Disks on Instability 

Real galaxies contain both gas and stars, and the system can be unstable even when 
both the stellar and gaseous components individually meet the requirements for stability 
(Jog & Solomon 1984; Elmegreen 1995; Wang & Silk 1994). Wang & Silk give an effective 
Q parameter for this two- fluid instability such that Qefi = «e//Q, where 

"e// = (1 + (/igfT,))"\ (12) 

In the limit of a low stellar surface density or large stellar velocity dispersion, we have 
«£// ~ 1-0; and Og// gets smaller as the stellar disk becomes more unstable relative to 
the gaseous disk. Following Wang & Silk, the parameters for the Milky Way disk yield 
= 0.72 in the solar neighborhood. 

How much should we expect this parameter to vary over our sample? The gas disk 
is a smaller fraction of the total gas mass in the earlier type galaxies, so the stellar disks 
are expected to have a more important influence on the stability of those disks. Scaling to 
the disk of the Milky Way (e.g. Kenney & Young 1989), the measured rotation velocities 
of these disks at Rhii imply average surface densities of 15 to 300 Mq pc~^. If the other 
parameters are held at their solar neighborhood values, then the highest surface density 
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disks should become unstable at a value of Ugas / fJ^crit that is ~ 25% lower than ag in the 
lower density disks. In other words, we would not be surprised to see values ranging 
from 0.25 to 1.0 across our sample. Better measurements of the stellar velocity dispersion 
and surface density would be needed to test whether the lack of such a trend posses a 
problem for the gravitational instability model. 



4.3. Consequences for Star Formation Prescriptions 

We have presented a picture where a single parameter, the gas surface density relative 
to the critical density for gravitational instability, determines whether the disk forms stars. 
Evolution in the cosmic star formation rate is likely driven, at least in part, by the evolution 
in the rate of galaxy - galaxy interactions with lookback time (e.g. Tan etal. 1999). In our 
star formation recipe, the environment influences the SFR indirectly through its impact on 
the large-scale gas distribution. Including the star formation threshold in the recipe allows 
the SFR to react non-linearly to small changes in the gas distribution. 

Our Ha atlas demonstrates that the distribution of star-forming regions in many disks 
is better described as asymmetric rather than axisymmetric. The correlations we find 
between local values of the instability parameter and the distribution of HII regions clearly 
indicate that the Toomre Q parameter provides useful guidance about the location of star 



forming regions in these disks. For example, in § |3.5.1| , we found that disks are more 
unstable on the side with the higher star formation rate. Prior to that we showed that star 
formation at sub-critical radii in NGC 2903, NGC 3031, and NGC 4548 is confined to the 
spiral arms. The average gas surface density was within a factor of three of the critical 
value in all three cases. The increase in gas surface density in the arms is expected to be 
larger than the increase in the epicyclic frequencey k (and hence critical density) caused 
by unresolved, streaming motions (Rand 1993; Elmegreen 1994), so the arm gas may be 
supercritical locally. Since these features will not be resolved in cosmological simulations 
either, it would be reasonable to lower the aq value in Eqn. 2 from 0.69 to ~ 0.3 in galaxies 
undergoing weak interactions to account for unresolved, super-critical structure in the gas 
distribution 

The fraction of gas residing in sub-critical regions of disks varies widely among the 
galaxies in our sample. Accounting for gas beyond the threshold radius produces the most 
significant correction to the global star formation rate. Among the spiral galaxies in our 
sample, the gas mass interior to the threshold radius Rhii varies from 0.14 to 0.94 of the 
total gas mass. Omitting the threshold radius from a model therefore leaves the global star 
formation rate uncertain by a factor of ~ 7. The median value of the gas fraction interior 
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to Rhii is 0.6. Ram pressure stripping has removed diffuse, atomic gas from at least these 
Virgo cluster galaxies - e.g. NGC 4394, NGC 4402, NGC 4569, and NGC 4689 (Warmels 
1988a, 1988b; Giovanelli & Haynes 1983; Kenney & Young 1989); yet a histogram of the gas 
mass fraction within Rhh for the Virgo subsample is indistinguishable from the histogram 
for the field subsample. 

5. Summciry 

Using a larger sample of spiral galaxies than previous studies, we tested the thesis 
that widespread star formation occurs where the gas disk is unstable to axisymmetric 
perturbations. We located the edge of the star-forming disk using our new Ha photometry, 
measured the gas surface density there, and derived the critical surface density for 
gravitational instability using the Toomre Q criterion. Our results confirm previous work 
by Kennicutt (1989) which found the threshold gas density varied by at least an order 
of magnitude among spiral galaxies but that the ratio of the gas density to the critical 
density was much more uniform. Among 26 galaxies with well-defined thresholds, we 
found the median ratio of gas surface density to critical density is ckq = 0.69 ± 0.2, with 
« (Rhh ) defined by eq. (2) and an assumed gas velocity dispersion of 6 km s~^ . At the 
thresholds, the gas is primarily atomic (molecular) in the disks with the lowest (highest) gas 
surface density. These results confirm the surprising accuracy with which a simple stability 
criterion describes the extent of the star-forming disk in large, objectively-defined samples 
of galaxies. 

This work also exposes the limitations of that simple model. The azimuthally-averaged, 
threshold values of the instability parameter, a (Rhii ), range from 0.3 to 1.2. A large part 
of that dispersion is due to uncertainties associated with asymmetric gas distributions, 
asymmetric mass distributions, and the extrapolation of the molecular gas surface density 
profiles. We demonstrated that the scatter in a (Rhii ) values is significantly reduced when 
Rhii and a (Rhii ) are fitted as functions of position angle. This result demonstrates that 
the simple Toomre Q stability criterion is useful for describing the distribution of star 
forming regions even in disks with quite non-axisymmetric gas distributions. The weak 
trend toward lower values of a (Rhii ) in Sab-Sbc galaxies, however, is not easily explained 
by any of the above systematics. These early-type spiral galaxies tend to have more massive 
stellar disks, so the lower values of aq would be consistent with the predicted de-stabilizing 
effect of a stellar disk on the gas disk. We find several sharp thresholds in the inner disks 
that are well described by gravitational stability alone. Star formation is also supressed in 
the centers of a few disks where the gas density is super-critical, but each of these disks 



-24- 



has a strong stellar bar. More detailed analysis of the star formation rate in the inner 
region of disks requires higher spatial resolution (cf. Jogcc 1999; Jogee etal. 2001). We 
conclude that streaming motions in the disk and the gravitational attracton of the stellar 
disk require small corrections to the simple instability parameter. Much of the uncertainty 
in current estimates of the instability parameter, however, should be resolved with complete 
two-dimensional mapping of the gas surface density and velocity field (e.g. Thornley et al. 
1999). 

One clear failure of the gravitational threshold model remains two sub-critical disks 
with widespread star formation in the inner disk - NGC 2403 and M 33. Both are low 
mass disks with low shear rates (i.e. slowly rising rotation curves). Elmegreen (1993) and 
HEB98 have suggested that star formation is hmited by the time available for clouds to 
grow via inelastic collisions (in the destructive presence of disk shear) rather than cloud 
formation, which they argue happens easily via magnetic instabilities. We demonstrate 
that their shear criterion offers an explanation for the high star formation rates in the inner 
disks of NGC 2403 and M 33 but find it does not consistently predict the outer (or inner) 
threshold radii. The vertical structure of these disks is also unusual in that the infall rate is 
high, and the high star formation rates may be indirectly related to this transitional state. 
Since the star formation rate depends more naturally on gas volume density rather than 
surface density, the impact of the vertical structure of gas disks on the star formation law 
remains an important issue in a more general sense. Unfortunately, radial variation in the 
gas velocity dispersion remains controversial due to the difficulty of the measurement (e.g. 
Ferguson etal. 1998; Sellwood & Balbus 1999). Among the galaxies in our sample, the 
intrinsic dispersion in the threshold value of the instability parameter is small for a constant 
velocity dispersion leaving little motivation for strong variation in velocity dispersion with 
radius. Given the difficulties in measuring velocity dispersions accurately, the issue could 
be settled by comparing the threshold radii reported here to models of the transition radius 
from disk self-gravity to halo self-gravity. 

The utility of the Toomre Q criterion extends beyond axisymmetric perturbations. 
Previous studies of sub-critical disks have focused on low surface brightness galaxies, which 
also have lower than average gas surface density. We showed that star-formation is similarly 
suppressed in disks with normal gas surface density if the rotation speed is fast enough to 
stabilize the disk. As the gas density approaches the critical density for instability, more 
isolated HII regions appear in the sub-critical disks. Similarly, in star-forming disks, HIT 
regions are found outside the threshold radius most frequently when the gas density in the 
outer disk is near (i.e. within a factor of three roughly) the critical density. Star formation 
is sometimes confined to a spiral density wave in these sub-critical environments, and the 
co-rotation radius of the spiral pattern may describe the outer edge of the star-formation 



-25- 



activity. This result seems to support the view advanced by Rudnick, Rix, & Kennicutt 
(2000) that tidal interactions impart a significant boost to the star formation rate in disk 
galaxies. 

Star formation thresholds in the outer disk change one's view of disk star formation in 
two ways. First, in a typical disk today, roughly 40% of the gas mass is not included in 
the empirical relation that defines the global Schmidt law for the SFR, and the analogous 
correction should be made when global star formation laws are applied in simulations. 
Second, disks are clearly not perfectly sclf-rcgulated if their outer disks are not forming 
stars. Significant velocity dispersions in these regions imply a physical process other than 
feedback from massive stars contributes to disk heating (Sellwood & Balbus 1999). Because 
of the threshold effect in the star formation law, the central density and angular momentum 
of the halo are clearly of fundamental importance to the star formation history of a disk 
galaxy (Dalcanton, Spergel, & Summers 1997; Mo et al. 1998; Mao et al. 1998) 
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Table 1: Data 



Galaxy 


Type'* 


d(Mpc) 




PA(°) 


Ha ^ 


Reference (d, RC, HI, CO) 


NGC 628 


SA(s)c 


11.4'' 


5° 


25 


K36 


1, 2 and 3, 1 


NGC 925 


SAB(rs)cd 


9.3 


55.2° 


102 


S90 


4, 2 and 5, 2, 1 


NGC 2403 


SAB(s)cd 


3.6 


60° 


125 


BS 


6, 2 and 5, 2, 1 


NGC 2841 


SA(r)b 


lO.O'' 


68° 


148 


K36 


1, 7, 7, 1 


NGC 2903 


SAB(rs)bc 


6.6'' 


60.0° 


44 


S90 


1, 2, 2, 1 


NGC 3031 


SA(s)ab 


3.6 


58.7° 


149 


BS 


8, 9, 10, 11 


NGC 4178 


Sdm 


17.0 


69.0° 


30 


S90 


12, 13, 14, 15 


NGC 4254 


SA(s)c 


17.0 


28° 


56 


S90 


12, 13, 14, 15 


NGC 4321 


SAB(s)bc 


16.1 


30° 


120 


S90 


16, 17, 14, 15 


NGC 4394 


(R)SB(r)b 


17.0 


25° 


108 


S90 


12, 18, 14, 15 


NGC 4402 


Sb 


17.0 


75.0° 


90 


S90 


12, 13, 14, 15 


NGC 4501 


SA(rs)b 


17.0 


58.0° 


140 


K36 


12, 13, 14, 15 


NGC 4535 


SAB(s)c 


16.3 


44° 





S90 


19, 13, 14, 15 


NGC 4548 


SBb(sr) 


16.1 


33° 


136 


S90 


6, 37, 14, 15 


NGC 4569 


SAB(rs)ab 


17.0 


65° 


19 


S90 


12, 13, 14, 15 


NGC 4571 


SA(r)cd 


17.0 


28° 


55 


S90 


12, 20, 14, 15 


NGC 4579 


SAB(rs)b 


17.0 


36° 


95 


S90 


12, 13, 14, 15 


NGC 4639 


Sbc 


17.0 


45° 


123 


S90 


12, 37, 14, 15 


NGC 4647 


SAB(rs)c 


17.0 


37° 


110 


S90 


12, 20, 14, 15 


NGC 4651 


SA(rs)c 


17.0 


42° 


71 


S90 


12, 18, 14, 15 


NGC 4654 


SAB(rs)cd 


17.0 


52° 


128 


S90 


12, 13, 14, 15 


NGC 4689 


SA(rs)bc 


17.0 


30° 


165 


S90 


12, 13, 14, 15 


NGC 4698 


SA(s)ab 


17.0 


57° 


164 


S90 


10 10 1/1 ir 

12. 18; 14, 15 


NGC 4713 


SAB(rs)d 


17.0 


49.0° 


100 


K36 


12, 20, 14, 15 


NGC 4736 


(R)SA(r)ab 


4.7'" 


35° 


122 


K36 


1, 21, 21, 1 


NGC 4826 


(R)SA(rs)ab 


5.4'' 


52° 


112 


S90 


1, 22, 22, 1 


NGC 5055 


SA(rs)bc 


8.4 


55.2° 


99 


K36 


1,2, 2,1 


NGC 5194 


SA(s)bc pec 


8.1 


20° 


170 


BS 


1, 23, 24, 25 


NGC 5236 


SAB(s)c 


4.r 


23.9° 


225 


BS 


26, 27, 27, 1 (36) 


NGC 5457 


SAB(rs)cd 


7.4 


18° 


44 


BS 


28, 29 and 30, 30, 29 


NGC 6946 


SAB(rs)cd 


4.8 


30.0° 


62 


BS 


21, 32, 33, 34 


NGC 7331 


SA(s)bc 


15.1 


74.9° 


167 


S90 


35, 7, 7, 1 



"Note. - Telescope used for Ha observation: BurrcU-Schmidt (BS); KPNO 36-in (K36); Steward 90-in (S90). 
''Published distance adjusted to Hq = 70km s"'^ Mpc"'^. 
^Distance to M83 group. 
''deVaucouleurs type from RCl 

^Inclination of the disk such that a face-on orientation is degrees. 

References. — (1) Young etal. 1995: (2) Wevers etal. 1986; (3) Kamphuis & Briggs 1992; (4) Silbermann 
etal. 1996; (5) Pisano, Wilcots, & Elmcgrccn 1998; (6) Graham etal. 1998; (7) Bosma 1981; (8) Freedman 
etal. 1994b; (9) Adler & Westpfahl 1996; (10) Rots 1975; (11) Sage & Westpfahl 1991; (12) Freedman etal. 
1994a; (13) Guhathakurta etal. 1988; (14) Warmels 1988c; (15) Kenney & Young 1988; (16) Ferrarese, L. 
etal. 1996, ApJ, 464, 568. (17) Sofue 1997; (18) Warmels 1988a; (19) Macri etal. 1998; (20) Warmels 1988b; 
(21) Bosma 1977; (22) Braun etal. 1994; (23) Sofue 1996; (24) Tilanus & Allen 1991; (25) Lord & Young 
1990; (26) Saha etal. 1995; (27) Tilanus & Allen 1993; (28) Kelson etal. 1996; (29) Kcnncy etal. 1991; (30) 
Bosma, Goss, & Allen 1981; (31) Sharina etal. 1997; (32) Carignan etal. 1990; (33) Tacconi & Young 1986; 
{?>A\ Ta.cconi Young- 1989: f35l Husrhes eta.l. 1998: fSfil Crosthwaite et a,l. 2000 nren^: Rubin etal. 
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FIGURE CAPTIONS 

Fig. 1. — Integrated Ha + [Nil] fluxes vs published values of Kennicutt & Kent (1983) and 
Kcnnicutt (1998). Open symbols denote galaxies with dependent calibrations - i.e. the data 
presented here were calibrated using the data plotted along the abscissa. 

Fig. 2. — (top) The Ha surface brightness profile for NGC 5236. The dashed line is the 
fitted exponential disk. The threshold radius is shown by the dotted line, (middle) Radial 
variation in gas suface density. Filled/open circles represent molecular gas density along the 
major axis from Young etal. (1995) and Crosthwaite & Turner (private communication). 
The squares show the azimuthally-averaged atomic hydrogen surface density (Tilanus & 
Allen 1993). (The H densities were multiplied by a factor of 1.4 to include the He mass.) 
The inner and outer spiral arms imprint the local HI maxima at 120" and 240" . The dashed 
line is the critical surface density for local gravitational instability, (bottom) Ratio of total 
gas surface density to critical density. 

Fig. 3. — Atlas of Ha + [Nil] images and radial surface brightness profiles. A vertical, 
dotted line marks the radius of the threshold in the azimuthally-averaged star formation 
activity. Dashed lines show variations in threshold radius with position angle. North and 
east are marked with bars 30" long. 
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Fig. 4. — Radial variation of the instability parameter in sub-critical disks. If no threshold 
radius was found, the radius of maximum instability was used to normalize the abscissa. 

Fig. 5. — Gas surface density in sub-critical disks. 

Fig. 6. — Mass fraction of atomic gas versus total gas mass at the star formation thresholds. 
The sub-critical disks are represented by open triangles. 

Fig. 7. — Azimuthal variation in the instability parameter. The solid lines represent an 
azimuthal average. The dashed lines are cuts along the major axis. Note that in (b) the HI 
surface density was measured at position angles of ±90°, and the molecular gas was measured 
along the major axis which is 60°off this HI strip. 

Fig. 8. — Radial variation in the instability parameter. The disks with the least azimuthal 
variation in threshold radius are shown. The galactocentric radius is normalized to the 
threshold radius. Dotted lines mark the mean ratio of the gas surface density to the critical 
density at the threshold radius. 

Fig. 9. — Radial variation in gas surface density. The disks with the least azimuthal variation 
in threshold radius are shown. 

Fig. 10. — Instability parameter at the star formation threshold vs. revised morphological 
type (RC2 T-type) of the galaxy. (Note that galaxy type is an integer quantity, and 
the decimals are used merely to separate points.) Filled (open) symbols represent surface 
densities dominated by atomic (molecular) gas, respectively. From left to right, the points 
represent these galaxies: NGC 4736 (2.0), NGC 4826 (2.1), NGC 4569 (2.3), NGC 4394 
(3.0), NGC 4402 (3.1), NGC 4501 (3.3), NGC 7331 (3.4), NGC 4579 (3.5), NGC 4689 

(4.0) , NGC 2903 (4.2), NGC 5055 (4.3), NGC 4321 (4.4), NGC 4651 (5.0), NGC 5236 

(5.1) , NGC 628 (5.2), NGC 4254 (5.3), NGC 4535 (5.4), NGC 4647 (5.5), NGC 2403 (6.0), 
NGC 6946 (6.1), NGC 4654 (6.3), NGC 925 (7.0), NGC 4713 (7.2), and NGC 4178 (8.0). 

Fig. 11. — Instability parameter normalized to the location of the threshold in the inner 
disk. The dotted lines shown the median value of the instability parameter at the outer 
thresholds for reference. 

Fig. 12. — Comparison of the Toomre Q and shear stabihty criteria in the inner disk. The 
observed gas surface density is shown by the solid line for comparison. The radial coordinate 
is normalized to the star formation threshold. 
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Fig. 13. — Comparison of the Toomre and shear instabihty parameters at the outer 
thresholds. The text describes why these six galaxies are shown. 
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